Analytical calculation of the solid angle 
defined by a cylindrical detector and a point 
cosine source with parallel axes* 
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Abstract 

We derive analytical expressions for the solid angle subtended by a right finite cir- 
cular cylinder at a point source with cosine angular distribution in the case where 
the source direction is parallel to the cylinder axis. As a subsidiary result, an ex- 
pression for the solid angle subtended by a disc detector at a spread disc source is 
also provided, in the case where the two discs have a common symmetry axis which 
is also coincident with the source direction. 

Key words: solid angle, point cosine source, cylindrical detector, cylinder, analytic 
expressions 
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1 Introduction 



In many situations in radiation physics the value of the solid angle sub- 
tended by a circular cylindrical detector at a point source is neede d. The 
case o f an isotropic p o int source has been treated to gr e at extent (I.Taffe 



(Tmsl. iManklinl (Kh . lMa.sket et aft dl95fih basket! dl^h. iGillesniej JUS 
Gardner and Verghesel (|l97ll h lOreen et aft f)l974l b IPratal (|2003bh b 



In a recent work ([PrataL |2003a) we derived analytical expressions for the solid 
angle subtended by a cylinder at a point cosine source, under the restriction 
that the source and cylinder axes are orthogonal to each other. In the present 
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work we obtain similar expressions in the case where the source direction is 
parallel to the cylinder axis. An expression (eq. 33) for the solid angle sub- 
tended by a circular disc with symmetry axis parallel to the source direction 
is also given as an auxiliary res ult. It should be m entioned that this latter 
expression appeared previously ( Hubbell et ak 1961 . eq. 29) in a slightly dif- 
ferent context. The authors detailed a quite general procedure to calculate 
the response of a small detector to an axially symmetric source with arbitrary 
polar angle distribution, uniformly spread on a circular disc. In that work, 
Q C i rc (here eq. 33) is interpreted as the response of a plane detector parallel to 
a Lambertian uniformly distribu ted disc source. In the s ame w ork the result 
is also credited to other authors ( Herman ( 1900t ). Footd ( 1915i l). 



With the on-going computer revolution, very complex problems otherwise im- 
possible to tackle, are routinely addressed using numerical methods. Calcu- 
lations such as the one presented here, amounting to a one-dimensional inte- 
gration, can be quickly specified and performed in a desk-top computer with 
minimal effort 2 . Nevertheless, we believe that in the few fortunate occasions 
where simple expressions exist in closed form, they deserve an interested look 
and are worth the effort of finding them. There are various reasons for this. 
No matter how careful the error-checking procedure followed, any numerical 
result should always be taken with some degree of prudence. In fact, except 
for the most simple situations, a great deal of effort is normally put into 
testing the result, rather than into obtaining it. Analytical results can be of 
some help here by providing reliable tests to computational algorithms. Also, 
a closed expression can sometimes be given an intuitive interpretation like 
the one provided for Q C i rc in section (2.3) and, at the very least, enables the 
usage of analytical tools to search for extreme values, regions of monotony, 
reliable approximations, asymptotic expressions, etc. Furthermore, the exis- 
tence of a closed expression for a point source reduces the complexity of the 
problem when considering the more realistic case of a planar spread source. 
The required four-dimensional integral can of course be reduced 3 to a double 
integral by using the exact expressions deduced here for a point source. In 
section (2.4) we illustrate this by obtaining an analytical result for the solid 
angle defined by a cosine source evenly distributed on a circular disc and a 
disc detector, in the very special situation where the two discs share the same 
symmetry axis which is also assumed to be coincident with the source di- 
rection. Reducing the dimensionality of a numerical integration is of obvious 



2 In fact, some of the results displayed in the present work were checked against a 
numerical integration. 

3 In the case of an uniformly spread source on a disc, a reduction to a double 
integral could also be achieved with resort to the exact expressions obtained by 



Hub bell et an (|1961T ) for the radiation field created by such source at an arbitrary 



point. The remaining double integral would then be performed over the detector 
surface rather than over the source surface. 
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advantage because the calculation is less prone to numerical errors and can 
be performed faster. In fact, when using quadrature formulae like Simpson's 
or Gaussian rules to calculate multidimensional integrals, the total number 
of points (Nt) scales like = N d , where N is the number of quadrature 
points used in each dimension and d is the dimensionality of the integral. The 
number of evaluations of the integrand can then become prohibitively high, 
regardless of the computing resources available, making necessary the resort 
to the Monte Carlo method, which, of course, has the drawback of a small 4 
conver gence rate (only as 1/ \fN T ) . This was discussed in a review by Jamesl 



(1980). who set (in 1980) the feasibility limit to the use of a ten-point Gaus- 



sian rule to five dimensions (with limited computer resources available) and 
to ten dimensions (with 'unlimited' computer resources). 

To illustrate the behavior of the solid angle sample plots are presented in 
section (3). 



2 Solid Angle Calculation 



The solid angle (ft sur f) subtended by a given surface at a point source can be 
defined as 



n surf = jj /(n)dfl, (l) 

directions 
hitting surface 

where f(Q)dQ is the source distribution. In the case of a point cosine the 
distribution is defined with respect to some direction axis specified by the unit 
vector k and it is given by /(fi) = (fi ■ k+ \ fl ■ k|)/(27r). The (27r) _1 factor 
ensures that < Q sur f < 1- Setting fi = f2 • k = cos(u), where v is the angle 
between the two axes, there results that f(ft) = {/V 71 " (p — 0)i (fi < 0)} so 
that the source only emits into the hemisphere around k. In the following we 
shall consider the situation of a right circular cylinder with axis parallel to k. 
The source is assumed to be at the origin of the coordinate system and the z 
axis is chosen both aligned with k and parallel to the cylinder axis. The solid 



4 It should be mentioned that there is always some value of d above which the 
Monte Carlo method is faster than any quadrature rule. This is so because the 
convergence rate of the quadrature rule is that of the one-dimensional rule. For 
instance, Simpson's rule error is proportional iV -4 or to N T ^^ d in d dimensions, 
which means Simpson's rule should, in general, be slower than Monte Carlo for 
d > 8. 
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angle is then given by 




(2) 



where 9 is the polar angle from the z axis; ip is the azimuthal angle in the xy 
plane and the limiting angles are to be determined from the conditions that 
\i = cos(#) > and that each included (9, </?) direction hits the surface. 

The solid angle Q of the whole cylinder can in general be decomposed according 
to fl = Q cy i + Q circ , where Q cy i and Vt circ are the contributions of the cylindrical 
surface and of one of the end circles. To calculate these quantities we refer to 
figs. 1 and 2, where, for simplicity, we assume that the cylinder and disc axes 
lie in the xz plane. To obtain Q cy i it is sufficient to consider the situation 
depicted in fig. 1, where one of the end discs is at the same plane as the 
source (i.e. z — 0). Let Q cy io(L,r,d) denote the solid angle in this case and 
tt circ (L, r, d) the solid angle defined by the disc. In the case of the disc we 
distinguish the situation shown if fig. 2, where r > d, from that where d > r. 
For economy this latter situation is also represented in fig. 1, with reference 
to the circle at the plane z = L. 




Fig. 1. Notation for Q C ylo and for £l c i rc (z = L) when d > r 



y 




Fig. 2. Notation for fi c j rc when d < r 



Using eq. 2 it follows that 
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(a) (b) 

Fig. 3. Definitions of (a) p± , 4>± and ip Q (d > r) and (b) p + , 4> + (d < r) 



< IML- r.(1) = 7T- 1 jf(l - pi(y?)/[L 2 + P 2 _(^)]) ^ 
o 

= 7r- 1 |L 2 /[L 2 + p 2 _(^)] dp = TrM-fo,) (3) 



Q circ (L, d>r) = n- 1 J (p 2 + (<p)/[L 2 + p 2 + (y?)] - pl(p)/[L 2 + p 2 (p)]) ^ 

o 

= ^-\J L 2 /[L 2 + p 2 _{p)} dp- J L 2 /[L 2 + p 2 + (p)] dtp) 



= 7i l [A_{y ) - A + {p> )\ 



(4) 



7T 

n circ (L,d < r) = tt' 1 J p 2 + (<p)/[L 2 +p 2 + (<p)} dip 

o 

7T 

= l-n- 1 jL 2 /[L 2 + pl(p)] dtp 



o 

-lr 



= 1-^^)] (5) 

where, from fig. 3, 

p Q = arcsin(r/d) (6) 



and 



P±( ( p) = d cos V 9 =t v r2 ~ (d sin y?) 2 . 
The required integrals A± can be expressed in terms of the integral 

r L 2 1 r 2 - d 2 
I{L, r, d, 0+) = / -(1 + ) , 

J L 2 + p4(0 + )2 p+(0+) 

where 

p + {4> + ) = \Jd 2 + r 2 + 2dr cos(0 + ) 
To proceed we begin by calculating 7. 

1 Calculation of I 



The integrand in the rhs of eq. 8 can be written as 

\ L 2 + <P-r 2 1 r 2 -d 2 1 

2^L 2 + rf 2 + r 2 l + mcos(0+) + d 2 + r 2 1 - n cos(0+) ^ ' 

where 

m = 2rrf/(L 2 + rf 2 + r 2 ) 

and 

n = 2rd/(d 2 + r 2 ) . 
The integration is straightforward, giving 



T L 2 + d 2 -r 2 1 . Il-m ,<j> + .. 

7 = L 2 + ^ 2 + r 2 Vl-m 2 TT^ tan( T )] 

+ ^ + r 2 Vl-n 2arCtanl /lT^ tan( l L)] • 
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Since 



(r 2 - d 2 )/(d 2 + r 2 ) = {VT^ri 2 (r > d); - V / T 3 ^ 2 (d > r)} 

and 



( L 2 + rf 2_ r 2 )/(L 2 + rf 2 + r2) = 

{l - m/ra(l + Vl-n 2 ) (r > d); 1 - m/n(l - y/l - n 2 ) (d > r)} (14) 
there results that 



1 — ml nil + y/l — n 2 ) , 1 — m , </>+,, 
I = ' ; ^arctan [J tan 

+ arctan[W ^j— ^ tan(^± )] ; (r > d) (15) 



and 



/ = 1 ~ mA f ~ arctanU/lI^tan^)] 
VI -m 2 [ Vl + m v 2 n 

— arctan[y y- j- — tan(-^)] ; (d > r) . (16) 



To express A + in the rhs of eqs. 4 and 5 in terms of /, a change of variable to <p + 
represented in fig. 3 is made, using 0+/2 = arctan [sin(</?)p + /( r — d + cos(</?)p + )], 
where p + = p+((p) is obtained from eq. 7. It follows that A + = / and, changing 
the integration limits, that 



A, (tt) = I\* (17) 

and 

A + (^ ) = J|- /2+ ^ . (18) 
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In the case of A-, the integration variable is first changed to 0_ shown in fig. 
3 and given by 0_/2 = arctan [sin(</?)p_/ (r + d — cos(</?)p_)], where, again, p_ 
is defined through eq. 7. Then, 



where 

p_(0_) = ^ + r 2_ 2drcos(0_) . (19) 
Making a further change to = n — 0_ yields 

n r T 2 1 r 2 — d 2 

where eqs. 19 and 9 were used to write 

p_(0_) = p_(7r -<f>) = V / cP + r 2 + 2drcos(0) = p + (0) . 

By comparison of eqs. 20 and 8 : 

A^ Q ) = -I\l /2+Vo . (21) 
Using eqs. 21, 18 and 17, eqs. 3 to 5 can be written as 

n cyl0 (L,r,d) = -n- 1 I\: /2+ ^ o , (22) 

Q circ (L, d>r) = -tt-VI W + I\l /2+iP °) = (23) 

and 

n circ (L,r > d) = 1 -tT 1 /^ . (24) 

Since J is discontinuous at d — r (see eqs. 15, 16), it should be clear that eqs. 
23 and 24 do not mean that Q circ (L, d > r) — Q C i rc (L, r > d) — 1. 

We now turn to eqs. 16 and 15 to evaluate the integrals. Since tan(0 + /2)|^ +=7r / 2+(po = 
[(d + r)/(d — r)] 1 / 2 = [(1 + n)/(l — n)] 1 / 4 and making use of arctan(^) + 
arctan(l/z) = 7r/2 , (z > 0), one obtains 
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n cy i (L,r, d) 



7r 1 (arctan 



1+n 
1 — n 



1 — m/n{\ — \J\ — n 2 ) 



arctan 



1 + m J 1 — n 
1 — m V 1+n 



(25) 



For Vt circ there results 



IWM > = 1/2(1 - '-"/f-f^) 



VP 



(26) 



MA d < r) = 1/2(1 - ^-m/^ + VT^?) ) 



(27) 



Defining f3 and 7 by 



tan(/3/2) = ^(l + n)/(l - n) 



(28) 



and 



cos(7) = m , 
eqs. 25, 26 and 27 can be recast as 



(29) 



tt cu , c (L, d>r) = 1/2(1 



1 + cos(/3) cos(7) . 
sin(7) 



(30) 



n /r 1 /on cos(ff) + cosfr) 

n circ (L,r > d) = 1/2(1 — - -r 77^ 

sm(7j cos(p) 



(31) 



and 



7T 



-I// 3 



Q cy io(L,r, d) 
1 + cos(/3) cos(7) 



sin (7) 



arctan 



cot(^)cot(|) 



) • (32) 



From eq. 28 it is seen that it/2 < (3 < ir and that, for d > r, (3 = ir/2 + (p . 
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As functions of L, d and r, eqs. 26 and 27 can be rewritten to give fW c in 
terms of a single expression: 

fi c „, c (L, d, r) = 1/2(1 - + ~ r ) , (33) 

(r 2 + d 2 + L 2 ) 2 - Ar 2 d 2 



which is equivalent to lHubbell et 3 jlQfill eq. 29), except for a factor of 7r. 



2.2 Special values and continuity 

First, the cases L = and d = r are discussed. From eqs. 11 and 12, it follows 
that 

L — > m — ► n , 
d-tr^{n-»lAm-» mi} , 
where mi = 2r 2 / (L 2 + 2r 2 ) . 
Then, using eq. 25, 

fi c?/ /o(£ -> 0, r < d) = Q cy io(m -> n) = 
Q cy io{L ^ 0, d -> r + ) = Q C2/W (n -> 1, m -> mi) = 1/2 
and ^cyzo is thus discontinuous when L — > 0, d — > r + . 
Starting from eqs. 26, 27, 



ft c „, c (L ^ 0, d -> r±) = n circ (n -> 1, m -> m x ) = 1/2(1 - ^f^) 

f2 C j rc (L — >■ 0) = {0 (d > r), 1/2 (d = r), 1 (d < r)}. One concludes that f2 circ 
is continuous whenever L ^ 0. 

Finally, for an infinite length cylinder at a skew position from the source (i.e. 
d > r) and with one end at the source plane there is the upper asymptotic 
value 



4 /l — n 

Q cy io{L -> oo,r < d) = f2 c ^ Q (m -> 0) = 1/2 - 2/7rarctan(*7 ) (34) 

y 1 ~\~ Tl 

= <Po/n ■ (35) 
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The last result can be checked directly using eq. 3 and can also be interpreted 
as the solid angle defined by an infinite length rectangle. 



2.3 Geometrical interpretation of Q C j, 



•circ 



Vt circ can be given a simple geometrical meaning by noticing that 



^circ = 1/2(1 - C0S(5)) , 



(36) 



where 5, shown in fig. 4, is the the aperture of the solid angle cone mea- 
sured in the plane containing the source axis and the disc center. Eq. 36 
is easily proved by writing 5 = arctan[(d + r)/L] — arctan[(rf — r)/L\ = 
arctan[2rL/(L 2 + d 2 — r 2 )], where in the last step we used arctan(zi) — 
arctan{z2) = arctan[(z\ — z 2 )/(l + Z\Z-if\. There results that cos(5) = (I? + 



d 2 - r 2 )/^{L 2 + d 2 - r 2 ) 2 + (2Lr) 2 . Since (L 2 + d 2 - r 2 ) 2 + (2Lr) 2 = (L 2 + 



d 2 + r 2 ) 2 — (2dr) 2 , it is seen that eq. 36 follows directly from eq. 33. 



2.4 The solid angle for a parallel coaxial disc source 

The solid angle defined by a spread cosine source evenly distributed on a 
disc and a coaxial, parallel, disc detector can also be obtained in closed form. 
Referring to fig. 5, we see that the solid angle is given by 



where R s and Rd are the source and detector radii; and Q C i rc (L, Rd, p)pd9dp 
is the contribution of an elemental area which can be obtained by putting 
r = Rd and d — p in eq. 33. Then, 



y 




Fig. 4. Definition of 5. 




(37) 




^(L 2 + P 2 -R 2 ) 2 + (2LR d ) 2 



-1 P=R t 



■S 



(38) 



Jp=0 
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A little algebra yields 



Q S (L, R d , R s ) 



R d 1 " 

R s 2m s . 




(39) 



or 



Q S (L, Rd, R s ) 



-R d tan(7 s /2) 
Rs 2 



(40) 



where 



sin(7 s ) =m s = 



2RdR s 



(41) 



V + Rl + Rl • 



y 




Fig. 5. Notation for the coaxial , parallel, circular source and detector set. 
3 Results and discussion 

In order to define the general position of the cylinder, let L 1: L 2 be the z 
coordinates of the end discs and set L\ = L 2 + length > L 2 . 

The solid angle Q(Li, L 2 ,r,d) can be calculated using eqs. 25, 26 27, 11 and 
12. The several situations to be considered are summarized in table 1. 

As examples, we consider two detectors of radius 1 and lengths 5 (L\ = L 2 + 5) 
and 10 = L 2 + 10). 

Plots of Vt for each detector as a function of L\ are shown in figs. 6 and 7, for 
different values of distance d. When d < r the curves for the two detectors 
are essentially the same except for a displacement equal to the difference of 
the detectors lengths. This is so because for Li smaller than the length of the 
detector [L 2 < 0), the source is inside the detector and Q — 1; for L 2 > the 
solid angle is that of the circle (VI = Vt circ ) and it has the same value for both 
equal-radius detectors. 

When d > r, Q increases with L\ as the detector is drawn from 'behind' 
the source. For the smaller distances (e.g. d — 1.5) Q quickly rises to the 
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Table 1 

Expressions for the solid angle of the whole detector. L±, L2 (L\ > L2) are the z 
coordinates of the end discs. 



Li, L2 


d, r 


n 


Li < 




0(«) 


Li > 0, L 2 < 


d < r 


1(6) 




d > r 


^cj//o(-^l,^ 


Li > 0, L 2 > 


d <r 


£lcirc(L 2 ,r, d) 




d> r 


Q cy w{Li,r, d) - n cy iQ(L 2 , r, d) + n circ (L 2 ,r, d) 



(a) the source only emits into z > 



(b) source inside the detector 

asymptotic value (see eq. 34) and remains almost constant until L\ is increased 
to values bigger than the length of the detector (L 2 > 0). Then f2 falls off as 
the detector is moved away from the source. It should be clear that Q(d > r) 
is the same for both detectors if L\ is smaller than the length of the shorter 
detector (i.e. L\ < 5). 




5 10 15 20 

Ll 



Fig. 6. Solid angle defined by a cylinder of radius 1 and length 5. 




Fig. 7. Solid angle defined by a cylinder of radius 1 and length 10. 
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